Derivation of continuum stochastic equations for discrete growth models 
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£NJ ■ We present a formalism to derive the stochastic differential equations (SDEs) for several solid- 

on-solid growth models. Our formalism begins with a mapping of the microscopic dynamics of 
, growth models onto the particle systems with reactions and diffusion. We then write the master 

Pj^ ■ equations for these corresponding particle systems and find the SDEs for the particle densities. 

Finally, by connecting the particle densities with the growth heights, we derive the SDEs for the 
height variables. Applying this formalism to discrete growth models, we find the Edwards- Wilkinson 
equation for the symmetric body-centered solid-on-solid (BCSOS) model, the Kardar-Parisi-Zhang 
equation for the asymmetric BCSOS model and the generalized restricted solid-on-solid (RSOS) 
• model, and the Villain-Lai-Das Sarma equation for the conserved RSOS model. In addition to the 

consistent forms of equations for growth models, we also obtain the coefficients associated with the 
SDEs. 
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In recent years, the study of nonequilibrium surface growth has attracted considerable interest in both analyti- 
cal and computational physics A number of discrete growth models and continuum stochastic equations have 
been proposed to describe the kinetic roughening properties of surface growth . By studying these models and 
equations, we classify them into universality classes according to their scaling behavior and associate the continuum 
I ' stochastic equations with the given discrete growth models. 

In general, two methods have been widely used to establish the correspondence between a continuum growth 
equation and a discrete growth model. One method is to use Monte Carlo simulations to obtain the scaling exponents 
from the discrete model and compare them with those of the corresponding continuum equation. The other is to derive 
the continuum equation analytically from a given discrete model. Computational methodologies have contributed 
significantly to our understanding of epitaxial growth over the past few years and continue to do sounabated. Analytic 
derivations include methods using the principle of symmetry Q or reparametrization invariance (9), and approaches 
starting from the master equation p0|-^3[. In particular, a systematic method proposed by Vvedensky et al. pi 
has been successfully applied to the derivation of the continuum growth equations directly from the growth rules of 
the discrete model for several solid-on-solid discrete models pj]-|l3|]. The derivation procedure of Vvedensky et al. 
, consists of two steps. First, the discrete stochastic equation is derived for the discrete growth model beginning with 
• the master-equation description of the microscopic dynamics of the discrete model. Second, the discrete equation 
is transformed into a continuous stochastic equation via regularization by expanding the nonanalytic quantities and 
replacing them with analytic quantities. In this regularization procedure, the step function is approximated by an 
analytic shifted hyperbolic tangent function, which is expanded in a Taylor series. As pointed out by Pfedota and 
Kotrla [O], the choice of regularization scheme for the step function is ambiguous. Thus, the coefficients in the derived 
continuum stochastic equation cannot be determined uniquely. 

In this Rapid Communication, we present a method for deriving the continuum stochastic equations from the 
discrete growth models. Our method can be applied to most of the models that are accessible via the method of 
Vvedensky et al. In addition to the derivation of stochastic equations consistent with the numerical solutions for the 
O ■ discrete models, our method predicts the coefficients in the stochastic equations. Our method begins with mapping 
of the discrete models onto reaction-diffusion systems with hard-core particles and sets up the master equation of the 
microscopic dynamics in the form of the Schrodinger equation. Next, borrowing the method introduced in Ref. p4| , 
we derive the corresponding Fokker-Planck equation, and then the stochastic differential equation is obtained. We 
apply our method to three discrete growth models: the body-centered solid-on-solid (BCSOS) model, the generalized 
restricted solid-on-solid (RSOS) model, and the conserved RSOS (CRSOS) model. 

The BCSOS model is one of the simplest microscopic growth models. Consider a surface built from square bricks 
rotated by 7r/4, without defects such as overhangs and vacancies. Each bond i contains a step 5, = ±1. The growth 
dynamics is as follows: Choose one of the columns at random. If this column is at the bottom of a local valley 
(Si-i = — 1 and Si = +1), a particle adsorbs with probability p (and nothing happens with probability 1—p). If it is 
at the top of a local hill (Si_i = +1 and Si = —1), a particle desorbs with probability q (and nothing happens with 
probability 1 — q). If it is part of a local slope (Si-\ = Si), nothing happens. This model can be mapped onto the 
problem of the asymmetric exclusion process (ASEP). The process describes particles that hop independently with 
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hard-core exclusion along a one-dimensional lattice with a bias that mimics an external driving force. By denoting 
the descending slope as the particle (A) and the ascending slope as the vacancy (0), the microscopic growth dynamics 
can be mapped to a stochastic dynamical rule of the ASEP: the only transitions allowed for the site with neighboring 
bonds (i — are 

A% — ► %A with probability p, 

%A -> A% with probability q. (1) 

This model has been studied extensively in the literature by Monte Carlo simulations |L5| and more recently it 
was realized that it can be solved exactly Jig] . Furthermore, Derrida and Mallick calculated the diffusion constant 
associated with fluctuations of the current in the limit of p ~ q and derived the corresponding continuum stochastic 
equation with the coefficients correct up to the lowest order Jl7[| . Krug et al. found the coefficients for the totally 
asymmetric, that is, q = 0, case j^]. Here, we apply our method to the ASEP to derive the corresponding stochastic 
equation. For simplicity, we set p — 1 and q = x (0 < x < 1). For x — 1 the system is symmetric, whereas for x — 
it reduces to the totally asymmetric case |f9j| . 

Introducing the annihilation and creation operators (fij's and a|'s, respectively) satisfying the mixed commutation 
relations explained in Ref. JlJ] and defining the state vector \ijj\t) = J2c t)\C), the master equation can be 
written as a Schrodinger-like equation: 

J^hM) = -%>;*>, (2) 

where P(C; t) is the probability for the system to be in a given microscopic configuration C at time t, and H , called 
the Hamiltonian, is an evolution operator expressed in terms of dj's and aj's. From now on, any operator will have 
a caret (e.g., a,b) and any symbol without a caret should not be confused with an operator. Occasionally, the same 
symbol is used to represent an operator and a density (e.g., hi and a,i). The Hamiltonian generating the time evolution 
of the BCSOS model is found to be 



(a 4 aj +1 + xa\a i+ ij . (3) 



Since the diagonal terms that make H stochastic have no role in our formalism, for simplicity, we omit these terms 
here and throughout this Rapid Communication. By involving the commutation relations between the Hamiltonian 
and some relevant operators such as ajaj and a\dia^a,j, and using the property of the projection state (-|(aj + a,i) = (-| 
PI , where (-| = ^2q{C\, we find the Kramers-Moyal coefficients |2(J corresponding to the above Hamiltonian: 

Ci = a t+ i + cii-i - 2cii - (1 - x) [aj+i(l - a,i) - a,i(l - a,_i)] , (4) 
= — [a,i + xa i+ i - (1 - xjdidi+i] (6i+i t j - Sij) 

+ [a,i-i + xai - (1 - x)a,ia,i-i] {8ij - 8i-ij). (5) 

Notice the absence of the caret on the a's. From the coefficients Ci and Cij, we write down the discrete stochastic 
equation 

+ & (6) 

where is the local density of the particle at i and (£i = CijS(t — t r ). Replacing <2j by (1 — V/i)/2 and taking 
the continuum limit, we obtain the continuum stochastic equation for the height variable h for the BCSOS model 
(and equivalently for the ASEP): 

dh 1 x 1 I 3y -1 , 1 . , . o , , . /_\ 

with (£(r; i)^(r'; t')) = 4(1 + x) j o(l — p)5(r — r')<5(t — t') (p is the stationary-state density, which is the same as the 
initial density). For the symmetric process (x = 1), the corresponding equation is the EW equation, whereas for the 
asymmetric process (x ^ 1) it is the KPZ equation. (Compare the coefficients with those in Refs. p7| , ^8[ .) 

Next we apply our method to the generalized RSOS growth model that was introduced by Neergaard and den 
Nijs pl| . They studied this model using an elegant mean-field type approach and derived the deterministic part of 
the KPZ equation. Our method is able to produce the stochastic (noise) part of the KPZ equation as well as the 
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same deterministic part as the method of Neergaard and den Nijs. The RSOS growth model describes the growth of 
simple cubic surfaces in which only monatomic steps are allowed. The heights at the nearest-neighbor columns can 
differ by only Ah = 0,±1 (the RSOS constraint). After choosing one of the columns at random, one particle can be 
deposited at the site i+ 1/2 with a probability 1 according to the height differences at i and i+ 1. Employing the 
same parameter as in Rcf . pl[ | , we map this model onto the two-species hard-core particle system with the following 
processes: 



00- 


AB 


Ph, 


00- 


> BA 


?«, 


AB 


-00 


Qh, 


BA - 


-00 


P«> 


0A- 


-> A0 


Ps, 


0B- 


■* B0 


q s , 


A%- 


-» 0A 


■ q s , 


B% - 


-> %B 


Ps, 



where a particle A (B) stands for the ascending (descending) bond and a vacuum represents the flat bond. The 
corresponding Hamiltonian is found to be 



H 



E 



Ph^M+i + qvb\al +1 + q h a,ib l+ i + p v bid, i+1 



+p s a\a i+ i + q s a. t a\ +1 + q s b\bi+i + p s hb\ +1 



(8) 



Following the same steps as above, we obtain the Kramers-Moyal coefficients C" and Cf? (a,/3 indicate either A or 
B) as follows : 



C 



A _Ph + q v 



Vi(v i+ i + Vi-i) + 



Ph - q v 



Vi(v i+ i - Vi-i) - q h aib i+ i - p v bi-\ai 
[a i+ i + a,i-i - 2a; + Oi(bi+i + bi-i) - &;(a i+ i + aj_i)] 
[(o,-+i - Oi_i)(l - 2a») - (a l+ i - ai-x)bi - a t {b. l+1 - h-i)] , 
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[vi(a i+ i + Oj_i) + Oi(«i+i + 
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where w,- 
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(9) 



(10) 
(11) 



1 — aj — 6; and Cj {CB ) is obtained by the exchange a «-> & followed by £ + fe «-> i — fc in C 4 (Q* ). 
By introducing the local slope D = a — b and the step density S = a + 6, we derive the deterministic parts of the 
stochastic equations for these two parameters: 



dD 
~dt 



= V 



f q (l - S) 2 + s d (l - S)S + \h g {S 2 - D 2 ) 



+ ^s m V 2 D + J(2 Sm - a s )V [DVS - SVD] , 

— = 2c.(l - S) 2 - -a s (S 2 - D 2 ) + ,s rf V [(1 - S)D] 

+ -h g [DWS - SWD] - c»(l - 5)V 2 S' 

+ -s TO V 2 5 - -a, [S'V^ - DX7 2 D] , 
2 4 L J 



(12) 



(13) 



where we have used the same notation for the parameters as in Ref. |2jJ ; c s = ph + q v , «s = Pv + Qh, s m = p s + q s , 
f g = ph — q v , h g = p v — qh, and Sd = p s — q s - The equation for the step density S is the same as that of Neergaard 
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and den Nijs. As pointed out in Ref. |2l|], the equation for the step density contains a mass term and the step density 

-l 



reaches its stationary value Sq = 1 + y / a s /(4c s ) after a characteristic time r s = (2y/a s c s ) . Thus, although 

there are two order parameters, only the local slope fluctuates at time scales larger than r s . At larger time scales, 
the step density does not behave as an independent dynamic variable. It follows local fluctuations in the slope of the 
surface: 



S = Sq + ^T s a s D 2 + t s 



VD + --. (14) 



Substituting Eq. (jlj) for S and identifying D = V/i, the equation for the height variable becomes 

dh „o, A 



.„ - . vV'h + ^tyhf + s, (15) 



where 



f g (l - So) 2 + «„(! - S )S + - A h g Sl 



]r s [2s d (\ - So) - h g S ][h g So - 4/ fl (l - S ) + 2s d (l - 2S )] + \s m (\ - S ) + ^ 



A = + ^r s a s [/i s S - 4/ s (l - S ) + 2s d (l - 2S )], (16) 

and (£(x,t)£(x / ,i / )> = £>^<S(:r - a/)5(* - *') with Z% = c s (l - S ) 2 + s ro S (l - So) + a s S 2 /4. This equation is the 
KPZ equation corresponding to the general RSOS model. To compare these coefficients with numerical work, let us 
consider the simple RSOS model introduced by Kim and Kosterlitz (KK model) B. The KK model corresponds to 
Ph = Pv = Ps — 1 and qh = q v = q s = 0. We obtain the corresponding coefficients Woo = 4/9, A = D{£ /(2u) = 2/3, 
A = —5/6. These values are consistent with the estimated values from a numerical study by Krug et al. fll8|] . 

Recently, a different growth model with a RSOS condition has been proposed and studied by Kim et al. 
Instead of rejecting the particle when the RSOS condition is not satisfied, this model allows the deposited particle 
to hop to the nearest site where the RSOS condition is satisfied. Thus, this model has the constraint of a conserved 
growth condition and is called the conserved RSOS model. The detailed derivation of Villain-Lai-Das Sarma (VLD) 
equation from the CRSOS model will be published elsewhere |^2[ . Here we only sketch the procedure and report the 
result. The procedure is similar to the one used to get the KPZ equation from the RSOS model except for some 
complicated calculation. After some algebra, we find the VLD equation 

^ = -0V 4 h + \V 2 (Vh) 2 + v , (17) 
at 

where v = (21-12v / 2)/2, A = (10-3\/2)/2, and (r}(x, t)r](x' , t')) = D m 6(x- x')S(t-t') with D m = (2\/2- l)/2. In 
deriving the above equation, we kept only the most relevant terms, and found that neither the EW nor the KPZ term 
exists in the growth equation. It is known, however, that higher order terms of the form V(Vft.) 2n+1 (n > 1) generate 
the EW term by the dynamic renormalization group |23j ]. Hence we should investigate the possibility of occurrence 
of these terms. Indeed, we found that the dangerous term of the form V(V/i) 2n+1 does not arise in the derivation of 
the VLD equation |2^]. Consequently, we concluded that the continuum equation of the CRSOS model is the VLD 
equation. 

Although the VLD equation was derived by Huang and Gu |l3| using the master-equation description with the 
regularization procedure proposed by Vvedensky et al., there were some ambiguities in choosing the regularization of 
the step function and thus the coefficients could not be predicted, whereas, in our derivation, we are able to predict the 
coefficients for the VLD equation corresponding to the CRSOS model. Unfortunately, however, there is no numerical 
method, to our knowledge, to find v and A for a microscopic model. Numerical studies up to now can argue only that 
A may be positive 0, which is consistent with our derivation. 

In deriving the VLD equation, we recognized the interesting aspect of the RSOS model. When we allow hopping 
processes only up to distance Iq, the stochastic equation for the height variable eventually becomes the KPZ one p2[ . 
This is contradictory to previously reported simulation results |24| . Kim and Yook studied the RSOS model with 
finite-distance hopping by Monte Carlo simulation and concluded that there is a phase transition at finite Iq from the 
KPZ class (Iq = 0) to the VLD class (Iq = oo). However, this conclusion seems to be a finite-size effect. In fact, we 
have confirmed our argument by carrying out a Monte Carlo simulation with sufficiently large system size for several 

k pi- 
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In summary, we have presented a formalism for deriving the continuum stochastic differential equations correspond- 
ing to discrete growth models. Applying the formalism to the BCSOS model, we derived the EW equation for the 
symmetric process and the KPZ equation with exact coefficients for the asymmetric process. The RSOS model was 
also studied with the general probabilities for possible processes. We derived the KPZ equation with a fluctuating 
noise part as well as the deterministic part, which is the same as the result of Neergaard and den Nijs. For the special 
case with p = 1 and q = (the KK model) , our coefficients are consistent with the numerical results of Krug et al. 
]18||. Finally, we applied our formalism to the conserved RSOS model. For the CRSOS model, we found that the 
VLD equation is the corresponding continuum SDE. However, if we allow only finite hopping (Iq < oo), the system 
belongs to the KPZ class eventually. We also predict the coefficients of the VLD equation for the CRSOS model. 
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